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Abstract 

The abundances of antiprotons and protons are considered within momentum-integrated Boltz- 
mann equations describing Little Bangs, i.e., fireballs created in relativistic heavy-ion collisions. 
Despite of a large antiproton annihilation cross section we find a small drop of the ratio of antipro- 
tons to protons from 170 MeV (chemical freeze-out temperature) until 100 MeV (kinetic freeze-out 
temperature) for CERN-SPS and BNL-RHIC energies thus corroborating the solution of the previ- 
ously exposed "antiproton puzzle" . In contrast, the Big Bang evolves so slowly that the antibaryons 
are kept for a long time in equilibrium resulting in an exceedingly small fraction. The adiabatic 
path of cosmic matter in the phase diagram of strongly interacting matter is mapped out. 
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I. INTRODUCTION 



The abundances of hadrons emerging from relativistic heavy-ion collisions can be de- 
scribed with surprisingly high accuracy by a thermo-statistical model Over a wide range 
of bombarding energies essentially two parameters, the temperature T chem and the baryo- 
chemical potential /x c /iem, adjusted to experimental data give a smooth curve T c / iem (/i c / iem ), 
the so-called chemical freeze-out curve, being an important landmark in the phase diagram 
of strongly interacting matter. At a given beam energy the variation of T c hem, Hchem for 
various centralities is fairly small Q]. 

The transverse momentum spectra of various hadron species, at a given beam energy, 
may also be described by a set of parameters T kin , jt% n , where additionally a flow parameter 



may be employed. It was argued, e.g. in [3j, and later confirmed, e.g. in that T c h em > Tkin 
holds in central collisions at high energies. This is often interpreted in a schematic picture 
as ceasing of inelastic (chemical) reactions at T c hemi while in a later stage of the fireball 
evolution a kinetic freeze-out happens at Tkm where the elastic interactions are no longer 
efficient enough to modify the momentum distributions. A different view was advocated 
in with Tchem ~ Tk in . In the former scenario there is still "hadronic life" after freezing 
out the chemical composition, while the latter scenario would mean a sudden termination 
of inelastic and elastic reactions roughly at the same time. 

A possible problem in the commonly accepted first scenario with Tchem > T^ n was pointed 
out in 6]: The antiproton (p) annihilation cross section is so large that the p survival 
until kinetic freeze-out needs special consideration. This was exposed as the "antiproton 
puzzle" . In it was shown that chemical off-equilibrium combined with multi-hadron 

reactions give rise to sufficient antiproton regeneration ensuring a moderate drop of the ratio 
of abundances p/p from T c hem until T^ n for CERN-SPS energies. This fits well in the above 
first hadro-chemical scenario 8]. 

It is the aim of this note to reconsider the off-equilibrium evolution of antiprotons in 
an expanding hadron fireball. We extend the consideration of 0] towards BNL-RHIC 
and CERN-LHC top energies. In a schematic model we show that the p/p ratio at CERN- 
SPS energies decreases from a given temperature T c hem ~ 170 MeV until kinetic freeze-out 
temperature Tj^ n « 100 — 120 MeV only by a small fraction, while at BLN-RHIC energies 
the variation of p/p in such a temperature interval is even negligible. We expect that at 
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CERN-LHC energies the latter statement also applies. We employ here a transparent (thus 
simplified) model to elaborate the features of the evolution towards off-equilibrium after 
an assumed chemical equilibrium which complies with the mentioned hadro-chemical model 
In contrast to this, the involved and sophisticated transport models give a much more 
complicated description, both for hadro-chemistry in general 9j and for the antiprotons 
especially [h], 
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Furthermore, we emphasize similarities and differences to antimatter (baryon) evolution 
in the early universe after confinement. The adiabatic path of cosmic baryon matter is 
exhibited in the phase diagram of strongly interacting matter. 

The paper is organized as follows. In Section 2 we discuss the appropriate evolution 
equations for protons (baryons) and antiprotons (antibaryons) coupled by a conservation 
law. The main emphasis is devoted to applications in relativistic heavy-ion collisions (Sec- 
tion 3) covering situations with baryons and antibaryons being nearly symmetric or very 
asymmetric. Section 4 describes the cosmic baryon matter. The summary can be found in 
Section 5. Appendix A sketches the derivation of the used evolution equations. 



II. OFF-EQUILIBRIUM EVOLUTION OF ANTIPROTONS 

Our starting point is the pair of momentum-integrated Boltzmann equations (see Ap- 
pendix A) 
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where Y± — n±/s are the baryon densities n± normalized to entropy density s = ^-h e ffT 3 . 
We attribute Y + ,n + to protons (p) and Y_,n_ to antiprotons (p). h e ff is the effective num- 
ber of degrees of freedom of hadrons in the fireball with comoving volume V(t). Here, we 
assume a spatially homogeneous fireball and apply for protons and antiprotons the Boltz- 
mann approximation, which is appropriate for the following since x = m^/T is larger than 
unity. Therefore, Y eq = -^^^x^K^ix) with the Bessel function K 2 . Both evolution equa- 
tions refer to pair-wise annihilations and regenerations, thus Y + — Y_ = rj = const, i.e. 
Y + = Y_. The two parameters in the evolution equations, and n, determine together 
with initial conditions the off-equilibrium dynamics. In what follows we are going to explore 



their interplay for Little Bangs and Big Bang. 

For Little Bangs (£ = 4, see Appendix A) one has 



A {4 ) = 3(av)fm 3 N ^h eff (l + r^j . (2.3) 

It encodes essentially the thermally averaged annihilation cross section (crv), and h e ff ac- 
counts for the time variation of the effective degrees of freedom. In deriving Eqs. (12.11 |2.2j) 
for Little Bangs we use adiabaticity of the expanding fireball, i.e. sV = const. Instead of a 
specific expansion model we utilize V/V = f with r as characteristic time scale. This allows 
us to formulate the evolution as a function of x instead of a function of time. 

In the Big Bang (£ = 2, see Appendix A and Q, the dimensionless factor A«) 



reads 



A (2) = (av) g]' 2 M Pl m N (2.4) 



with gV 2 = h e jf g e jf 2 f 1 + ^ J an< ^ determining the effective degrees of freedom 

h9eff 



relevant for the energy density, i.e. e = ^rg e f/T 4 , and M P i as Planck mass and uin again 



as nucleon mass. Clearly, different h e ff's apply in Little Bangs and Big Bang. 

Using (av) = C/m 2 with C = 0(1) and with m n denoting the pion (tt°) mass one arrives 
for f ~ 5 fm/c at A( 4 ) ~ (9(10 4 ) as an estimate (for h e ff see below), while Am ~ C(10 21 ) 

1 /2 I 

for gj ~ 4 [16[ highlights the vast difference of Big Bang and Little Bang dynamics. 



Actually, the thermally averaged cross section is (av) = ^ 2x ^ ^j.^a^ 1 ^^^'"^ 13] with 



= T£Ve-4a;7(2x), see also (nj. Employing a(p lab ) = (4Qp^ v/c + 2Ap^ GeV/c ) mb 
[2] or (38 + ^PiabGeV/c) f° r the annihilation cross section one arrives at (av) ~ 42.5 
- 47 mb in the temperature interval 170 - 100 MeV or 51.5 mb fairly independent of T, thus 
yielding C ~ 2. 



III. RESULTS FOR LITTLE BANGS 

Equations (12. 1[ 12.21) are of Riccati type with no general analytical solution in closed 
form. The solutions, to be found numerically, depend on A^ and the initial condi- 
tions (encoded in rj). For the latter ones we employ T chem = T = 170 MeV and 
fJ'chem = yUo = 250 MeV (SPS top energy resulting in n,_/n+ = 0.052) or fi c hem = 
/i = 25 MeV (RHIC top energy resulting in n_/n + = 0.75), thus neglecting a possi- 
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ble small change of chemical freeze-out temperature when going from SPS to RHIC but 
catching the typical values extracted in the data analysis Q, SPS energies are in 
the realm fi c h em > T chem , while RHIC and LHC operate in the region /i c /i em < T chem . 
A good approximation of h e ff is provided by h e ff(T) — hi + b(T — 1\) + c(T — 
T x f with b=[h-h 1 - (h 3 - h)(T 2 - T x y/{T Z - T x f\ / [T 2 -T x - (T 2 - T X ) 2 /(T 3 - T x )\ 
and c = (h 3 - h x )/(T 3 - T x ) 2 - b/(T 3 - T x ) and (T X ,T 2 ,T 3 , h x , h 2 , h 3 ) = 
(0.100, 0.150, 0.175, 3.0, 9.8, 17.5) for small chemical potentials as appropriate for RHIC top 
energies and (0.100,0.150,0.175,5.411,13.150,21.750) for an isentrope with entropy per 
baryon equal to 23 as appropriate for SPS top energy. These values are for the resonance 
gas model with the first hundred hadronic states as used in [l| and for temperatures in GeV. 
While the scales (av) and f as well as h e ff are condensed into one dimensionless parameter 
A(4) which would be useful for h e ff = const, the strong variation of h e ff in the consid- 
ered temperature range along isentropic trajectories makes A( 4 ) not a concise characteristic 
quantity. Therefore, we present the results for various values of C to expose the influence 
of annihilations and to discriminate them from diminishing densities due to expansion. 

In Fig. [T]the change of n^/n + (denoted by p/p) as a function of the temperature is ex- 
hibited. We note that, similar to the consideration in [7], the ratio drops only by a small 
fraction even when extending the evolution down to a temperature of 100 MeV. The reason 
is the small change of Y±. The experimental mid-rapidity value of p/p is 0.058 ± 0.005 (sta- 
tistical error) in central collisions of Pb + Pb at E^ eam = 158 AGeV [18| . Clearly, our above 
mentioned initial value is already somewhat below this experimental value. If one would 
change [iq to 235 MeV, being still in the range of admissible values according to [1], the 
start value would be n_/n + = 0.063 and all curves exhibited in Fig. [1] are then upshifted, 
roughly by factor 1.2. As a consequence, values of C < 4 would be compatible with the 
data. While this range of C is realistic, such a fine tuning is not appropriate given the 
schematic character of our model (use of a characteristic time scale f for the expanding 
fireball, no feeding etc.). This refrains us from considering further details like the centrality 
dependence. Nevertheless, one may consider the beam energy systematics for central colli- 
sions. Instead of individually selected values of fi c hem and T chem , one may use a global fit of 
many hadron ratios with the parametrization [1], ^i c hem — 1303 MeV / (1 + 0.286y / J/v r /v/GeV) 
andT chem = 162 MeV (l - (0.7 + exp[( v /^7/GeV - 2.9)/1.5])- 1 ) (where is in GeV) 

and get initial values for n_/n + being 15%, 1.5%, 18%, 30% and 29% larger than the ex- 
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perimental ratios p/p quoted in [18j| for central Pb + Pb collisions at SPS beam energies of 
158, 80, 40, 30 and 20 AGeV. The offset of ri-/n + above the experimental p/p value [3] 
at ^/snn = 200 GeV is 36% according to this parametrization. With the exception of the 
data situation for beam energy of 80 AGeV, there is room for a 20% drop, in average, of 
p/p towards kinetic freeze-out, consistent with the results in Fig. [TJ 

In contrast to the ratio p/p, the densities rapidly change with dropping temperature, 
as exhibited in Fig. [2j A notable point is the increasingly strong departure of n_ from 
the chemical equilibrium value n eq towards kinetic freeze-out. For this comparison we have 
determined n± by sY eq exp(±/x/T) with p from the resonance gas model along the isentropic 
curve imposing baryon conservation (a convenient parametrization is provided by p = a/T x + 
b with (x, a, b) = (1, 0.045198563, -0.015873836) - T and p in units of GeV). The emerging 
relation n eq > n + (see Fig. [2]) is surprising at the first glance. In fact, baryon conservation 
enforces > npV{t = 0)/V(t > 0) for the net nucleon density. Since the antibaryon 
density in the situation at hand is significantly smaller than the baryon density, the nucleon 
density also decreases slower than V(t = 0)/V(t > 0). As our evolution equations impose 
pair-wise annihilations and regenerations, for n + 3> n_, n+ goes approximately with V(t = 
0)/V(t), while the fiducial density obeys n eq > n+V(t = 0)/V(t > 0). This ostensi 
ambiguity can be resolved by introducing effective chemical potentials in line with 
This, however, is not necessary for the present purposes, as the combination n eq n eq (or 
Y+ q Yl q ) enter the evolution equations, and neither n eq nor n eq separately: In the employed 
Boltzmann approximation, the chemical potentials cancel in Y+ q Yl q = Y^ q . Discarding 
this gain (or recombination) term would result in a significantly reduced ratio p/p, thus 
emphasizing the importance of the back-reaction, as already stressed in 

m. 

(In |7| the 

regeneration term is of utmost importance to counteract the stronger annihilations at the 
total baryon content.) 

Summarizing, the inspection of Fig. [2] reveals the strong deviation from equilibrium, i.e., 
despite of the large annihilation cross section the expansion is too rapid to maintain chemical 
equilibrium of antiprotons. For larger values of {av) parameterized by C, n_ follows more 
closely the fiducial density n eq and, as a consequence, the antiproton to proton ratio drops 
stronger during cooling (see also Fig.[T]). A larger value of the expansion time scale f reduces 
also the ratio p/p: For f = 10 fm/c we get 0.040 at T = 100 MeV. In contrast, a shorter 
expansion time scale keeps the ratio at higher values, say 0.048 for f = 3 fm/c. As n + is 
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essentially not modified, one can infer the corresponding values of n_ from these numbers. 

Reference l| (first quotation) states that the thermo-statistical model applies not only 
to selected ratios of hadron yields but also to yields themselves. It is instructive, therefore, 
to consider the evolution of the yields from chemical freeze-out towards kinetic freeze-out 
temperatures. We normalize the yields according to (niV)i 00 Mev/( n iV)i70MeV an d exhibit 
in Fig. [3]the dependence on C = (av)/m%. The yield of antiprotons depends smoothly on C; 
values C > 3 would cause a sizeable annihilation which would obscure the hadro-chemistry 
picture. 

For RHIC energy the same features hold. However, different initial conditions cause some 
different evolution: The ratio p/p is fairly insensitive to C, see Fig. HJ The departure from 
equilibrium for protons is nearly as strong as for antiprotons, as shown in Fig. (Here, the 
above parametrization of /x(T) with (x,a,b) = (4.8,2.2671929 x 10" 6 , 0.013797031) applies.) 
Also the yields in Fig. [6] drop with increasing values of C as they are more dragged by 
the respective equilibrium values. For C < 10 the reduction via annihilation is less than 
10% thus not invalidating the consistency of hadro-chemistry with a late kinetic freeze-out. 
Experimentally, one finds p/p = 0.731 ± 0.011 ± 0.062 19] in central collisions Au + Au at 
y / J/v r /v = 200 GeV which is fairly independent of centrality. This value compares well with 
the results in Fig. HI even for large C. 

For smaller values of \x c hem the difference of n_ and n + becomes smaller: Annihilation 
diminishes both n_ and n + by the same amount thus keeping the ratio n_/n + nearly con- 
stant. This consideration applies in particular for LHC In contrast, for the above baryon- 
antibaryon asymmetric situation at SPS, given by a larger value of fi c hem, the annihilation 
of p by the small admixture of p is not severe, and only the evolution equation for ri_ is 
sufficient, as exploited in [i]]. 

Resolving the "antiproton puzzle" means demonstrating that the p/p ratio does not 
change noticeably from T c h em until T^ n . At high beam energies (say, for RHIC and LHC 
energies) corresponding to smaller values of the normalized particle-antiparticle number 
difference 77 this seems to be fairly robust. At SPS one sees already a sensitivity to the 
interplay of thermally averaged annihilation cross section and expansion dynamics. At AGS 
and later on planned FAIR energies one expects a stronger drop of the p/p ratio. However, 
the differences of T c h em and Tki n may be much smaller so that again annihilation is less 
severe. 
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Our results base on a few assumptions which we recollect here: (i) kinetic equilibrium 
(which is left at the kinetic freeze-out, see for dealing with the freeze-out itself); (ii) 
spatial homogeneity; (iii) expansion dynamics characterized by one time scale; (iv) restriction 
to one hadron species, p,p (the other hadrons are implicitly in the heat bath, encoded in h e ff 



351]); (v) chemical equilibrium at T chem , ii c hem\ ( y i) detailed balance and unitarity. Items (i - 
v) are relaxed in transport codes, which also attempt to include (vi). We insist here to arrive 
at a qualitative and transparent understanding of an aspect of chemical freeze-out. Items (i) 
and (ii) are related with the derivation of the employed form of the momentum integrated 
Boltzmann equation, see 3] for details. Item (iii) refers to the fact that we consider a 
specific expansion pattern of the fireball (which is assumed to contain a homogeneous matter 
distribution according to item (ii)) characterized by V/V = f with f as expansion time scale, 
see Appendix A. Item (v) fixes the initial conditions in agreement with the thermo-statistical 

n 

model pj. 

Strictly speaking, item (vi) refers to the balance p + p^X + Xas the underlying 
Boltzmann equation including a binary collision kernel leads to our evolution equations 
(12.11 12. 2p . One may think, however, that both X and X represent clusters of pions. Such 



states X, X are successfully considered in 



21] as two-meson doorway states coupling in 



turn to multi-pion states. This approach, also known as minimal two-body model, realizes 
the nearest threshold dominance and describes the pp annihilation data in the energy region 
relevant for our purposes. It bridges to |f3,Q], where such multi-pion collisions were considered 
as key to resolve the antiproton puzzle together with chemical off-equilibrium in the pion 
component. As further possible ingredients one may consider the role of baryon excitations 
(as sources for further annihilations and regenerations) and feeding of the ground state 



baryons eventually observed, as done in 



22| for RHIC energy. 



IV. ANTIBARYONS IN BIG BANG 

Due to the numerically large value of A( 2 ), the evolution of n± follows closely the equi- 
librium values n± for a long time. In other words, Y± are dragged by Y± q as evidenced 
by solving numerically Eqs. (2.1, 2.2). In the temperature region T > 1 MeV, Y± « Y± q 
represent highly accurate solutions of the evolution equations. Freeze-out of the antinucleon 



annihilation happens at temperatures being considerably lower than the MeV scale 15) thus 
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yielding an exceedingly small antibaryon density in the assumed homogeneous scenario (for 
inhomogeneous scenarios cf. j^)- Indeed, using n± m n± = sY eq exp{±/i/T}, one finds 
a rapid dropping of the scaled antibaryon density, represented by n_/T 3 , with decreasing 
temperature, see top panel of Fig. [7J (/j, is determined by Eq. (4.1) below). For the given 
small surplus of baryons, the scaled density n + /T 3 follows closely n_/T 3 until 40 MeV; on 
the exhibited scale in Fig. [7J the difference of n + /T 3 to n_/T 3 is not visible. Below that 
temperature of 40 MeV, however, n + /T 3 stays approximately constant at ns/T 3 ; the still 
continuing annihilations diminish n + /T 3 only marginally 36j since n_/T 3 became small. 
The turn of n + /T 3 from dropping to the approximately constant value near ns/T 3 at tem- 
peratures of 40 MeV is entirely determined by the value of r\ 24| which quantifies the baryon 
surplus, encoding also the chemical potential. 

In the cosmic evolution, say after confinement, rj ~ 10~ 10 is presumably realized, governed 
by the observed ratio of baryons to photons being 6.12 x 10~ 10 (ACDM 3- year WMAP-only 



data 



25]) and relying on adiabaticity. This ratio also determines the adiabatic path of 



cosmic matter after confinement until the onset of primordial nucleo-synthesis at T ^ 1 
MeV, see Fig. E] which is based on baryon conservation expressed by 

Strictly speaking, this equation applies for h e ff = const and for situations where nucleons are 
essentially the carriers of baryon charge. Accordingly, the cosmic baryo-chemical potential 
fi evolves from 10~ 6 MeV at T = 170 MeV towards ~ 938 MeV prior to nucleo- 



synthesis at T < 1 MeV 



24 



261 ] . Features of Eq. (14. ip are: (i) min(fi) = 0(171^1]), (ii) 



max{n) = rriN, (in) crossing the /i = T line at T w itln/{— log^+- • ■) (this is the point where 
the difference of n + and n_ becomes large with dropping temperature: n + /T 3 then stays 
approximately constant at r](27c 2 /A5)h e ff, while ra_/T 3 continues dropping exponentially, 
see Fig. [7J and 0]), (iv) before the region /i ~ T is reached, the temperature drops as 
T oc m^/(log(fi/ 171^7]) + ■•■). The • ■ ■ in items (iii) and (iv) are for subleading terms. The 
variation of h e ff with temperature, according to the resonance gas model, pulls down the 
isentropic curve (solid line) for ^ < 1 keV so that the crossing with the dashed line occurs at 
a temperature of 165 MeV (instead of 225 MeV). For fi > 1 keV, the results of the resonance 
gas model with adiabaticity and baryon conservation are on top of the solid line. 

For an orientation, in Fig. [8] also chemical freeze-out points in relativistic heavy-ion colli- 
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sions from the analysis in [l| (first quotation) are exhibited. At the freeze-out temperature 
of about 160 MeV in Little Bangs at RHIC the baryo-chemical potential is about 10 7 times 
larger than in Big Bang. Otherwise, when the baryo-chemical potential values of Little 
Bangs at RHIC is achieved, the temperature in Big Bang is about 40 MeV (see also {27)]). 
being surprisingly high. 

During the evolution of matter after electro- weak symmetry breaking at T ~ (9(100 GeV) 
down to confinement at T c ~ (9(200 MeV), the strongly interacting matter dominates by 
far the pressure, the energy density and the entropy density [28]. The masses of carriers of 
baryon charge change in the confinement transition. Below 160 MeV the energy density is 
dominated by electro- weak matter for a long time including primordial nucleo-synthesis, see 
bottom panel of Fig. [7J Pions are exceptional, as they are sizeable in number and energy 
density contribution down to 5 MeV. The bottom panel of Fig. [7J evidences also the slow 
relative increase of the scaled baryonic energy density; much later (before "recombination") 
it will exceed the electro-weak matter thus turning the radiation universe into a matter 
dominated universe. 



V. SUMMARY 



In summary we contrast the baryonic antimatter evolution in Little Bangs (i.e. bulk mat- 
ter of fireballs created in relativistic heavy-ion collisions) and Big Bang from confinement 
towards the onset of primordial nucleo-synthesis. To expose the similarities and differences 
we focus on antiprotons in heavy-ion collisions and nucleons in the cosmic evolution within 
a schematic kinetic description neglecting the explicit coupling to other hadron states. The 
vast differences of Little Bangs and Big Bang are impressively described by the huge differ- 
ence of the dimensionless parameters A(£) governing the chemical freeze-out dynamics. In 
addition, the scaled net baryon densities encoded in r\ are also drastically different. 

In Little Bangs, the ratio p/p drops only by a tiny amount from the established chemical 
freeze-out temperature (which is defined by a multitude of other hadron abundances) until 
kinetic freeze-out for RHIC conditions, where antiprotons appear in a sizeable fraction. For 
a larger asymmetry of antiprotons to protons (as realized for SPS conditions) the ratio 
p/p changes still by a sufficiently small amount to maintain the consistency of the thermo- 
statistical model for hadron chemistry. The thermo-statistical model is an important tool 
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for mapping out the phase diagram of strongly interacting matter in the T — fi plane. In 
so far, it is important that its features conform with a detailed description of heavy-ion 
collisions. 

In homogeneous scenarios for the Big Bang, the chemical equilibrium is maintained for 
a long time thus resulting in an exceedingly small fraction of primordial antibaryons jl5 |. 
An interesting point is the adiabatic path of cosmic baryons through the phase diagram of 
strongly interacting matter and the related temperature of about 40 MeV (determined by 
if) during the orders-of-magnitude change of the chemical potential turning cosmic baryon 
matter from fi/T -C 1 into fi/T 3> 1. 
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APPENDIX A: MOMENTUM-INTEGRATED BOLTZMANN EQUATION 

Let be LjiVj = Cj the Boltzmann transport equation for particles of species i, where Lj 
is the covariant Liouville operator, Li = Pi-£^ — r^pfpf ^ with affine connection 
and Greek indices running from 0-3 (Einstein's sum convention applies only for them). 
Cj may be a general source term modifying the free-stream of the distribution function 
Ni(x,p) by collisions and fee ding by decays. A consequence of the structure of the Liouville 



operator is iV" = / LjA^IIj [29|, |30j and thus with the Boltzmann equation one gets iV" 



J Cjllj, where the coordinate- independent momentum element for particle of species i is in 
the notation of [30J IL = di(2Tr)~ 3 ^—g^^ with the fundamental determinant —g, to be 
built of the metric tensor g^ u , and particle degeneracy df, the semicolon stands for the 
coordinate-covariant derivative. The particle momenta p% are normalized to rest masses rrii 
as g^uPiPi = — m "i- The particle current is defined by N? = J Nip^Hi. For an observer 
moving with four- velocity v a the corresponding particle density is n« = —VaN®. 

We consider here bulk matter in sufficiently small volume elements where matter looks 
isotropically and homogeneously. Moreover, the matter is provided to be in collisional (i.e., 
thermal) equilibrium so that a common flow field u a (x) can be attributed to all particle 
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species. This means, N" = riiU a 3l]. Synchronizing the observer velocity with the flow by 
introducing a comoving coordinate system where u a = g$ = 8$, the l.h.s. of the balance 
equation for the density n i; (nj-u a ) ;Q = J CjlL, becomes ^= (a/ - S'^fl'o ),« = 77=7i(V~9 n ^' 



where the dot means time derivative in the respective comoving coordinate system. Intro- 
ducing with the reasoning of [30[] a comoving three- volume V, defined, e.g., by (n con V)' = 
for a conserved charge density n con or by entropy conservation, (sV)' = 0, the l.h.s. of the 
considered balance equation may be written as y(riiV)'. 

Binary elastic collisions, % + i' «-> i" + i'", cancel out the r.h.s. of the balance equation as 
they represent collisional invariants. For the moment being, we focus on binary annihilation 
processes, % + % <-> X + X [32j. In the absence of Bose condensation or Fermi degeneracy 
effects and with the assumption that species X and X have equilibrium distributions, the 
collision term / Cjllj can be evaluated (see 



13. 



32] for details) to result in the master type 



equation 



rii + rii 



V 
V 



-{av)(n i n- i - nfnf^ 



(Al) 



and an analog equation for i by the replacement i — > i. For Robertson- Walker- Friedmann 
cosmology, where strict isotropy and homogeneity is applied, this is an often derived and 



applied equation, see [13| 



Him 



32] , and is referred to as momentum- integrated Boltzmann 



equation. In heavy-ion collisions such an equation, supplemented by decay terms, is used too, 
for instance in 22|. In order to apply the above arguments, however, special flow symmetries 
must be required. We consider examples below. (Other approaches [33] to rate equations 
use a spatial average over the expanding fireball, where a Lorentz factor for transforming 
the time coordinate to the external coordinate system occurs additionally.) 



We follow further the arguments in 



13 



32] and reformulate the rate equation (1A1I) . 



instead as time evolution equation with respect to comoving or proper observer time, as 
evolution equation with respect to temperature T: Define the yields Yi = rii/s, use the 
definition of the comoving volume V by (sV)' = 0, and introduce the new variable x = mjq/T 
to arrive at 



heff T 



(A2) 



dx N ' V 1 1 * J x 4 V 45 JJ \ h eff 3TJ 
Clearly, further evolution equations are still needed describing the expansion dynamics of 
the system under consideration. For our schematic discussion in heavy-ion collisions we 
put V/V = f with f as characteristic expansion time. If the flow pattern is assumed to be 
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specific, the relevant time scales are henceforth prescribed. Examples for simple flow patterns 
are the four-velocities (i) = (chf), 0, 0, shr)), (ii) uF 1 = (chf), shr) cos0, shr) sin 0, 0), (hi) 

= (chf), shr) sin 9 cos0, shr} sin 9 sin 0, shr) cos 6>) with r) = | log((t+C) / (t— C))> f° r (i) C — 
(ii) £ = y/x 2 + y 2 and (iii) £ = \Jx 2 + y 2 + z 2 describing (i) purely longitudinal (Bjorken), 



(ii) axial-symmetric transverse und (iii) spherical expansion yielding with r — yi 2 — ( 2 for 
(l/V)(dV/dr) (i) 1/r, (ii) 2/r, and (iii) 3/r (here, i,x,y,z are Cartesian coordinates in 
Minkowski space-time where u a refers to; 9 and are usual cylinder or polar coordinates). 
The approximation leading to Eqs. (2.1, 2.2) consists in the replacement of the latter ratios by 
the inverse of a characteristic time scale, 1/r. Executing the transformation to the comoving 
coordinates via u a = dx d ^ yields u a = g^; the line elements in the comoving coordinates 
read (i) ds 2 = -dr 2 + dx 2 + dy 2 + T 2 df] 2 , (ii) ds 2 = -dr 2 + r 2 dfj 2 + T 2 sh 2 fjd(j) 2 + dz 2 , (iii) 
ds 2 = —dr 2 + r 2 dr] 2 + r 2 sh 2 r]d9 2 + r 2 sh 2 r] sin 2 9 d<p 2 representing time-orthogonal (Gaussian) 
coordinates, as required to arrive at Eq. (lAlj) . 

In Robertson- Walker- Friedmann cosmology with the line element ds 2 = —dt 2 + R 2 (t)dx 2 
the comoving velocity is as mentioned above; R(t) is the scale factor. The Einstein equations 
govern the dynamics. For conformal flat three-space and without cosmological constant one 
has 



3j V/V = 3R/R = 3JfG N e with Newton's constant G N = Mpf and the energy 



density e = ^^//T 4 . Combining these quantities appropriately one gets Eqs. (2.1, 2.2) 
with (2.4) from (IA2I) as compact notation of the evolution equations for protons or nucleons 
(i — > +) and antiprotons or antinucleons («—►—). 
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FIG. 1: (Color online) The ratio of antiprotons to protons as a function of temperature for various 
values of C = 1, 2, 10 (from top to bottom). Initial conditions for SPS as described in the text. 
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FIG. 2: (Color online) The proton density and antiproton density as a function of temperature. 
Initial conditions for SPS. Solid curves: actual densities n± for the off-equilibrium evolution (i.e. 
solutions of Eqs. (|2.1l I2.2|) ). dashed curves: fiducial equilibrium densities n± which need special 
explanation (see text). For C = 2. 
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FIG. 3: (Color online) The proton yield (upper curve) and antiproton yield (lower curve) at 100 
MeV normalized to those at 170 MeV as a function of C. Initial conditions for SPS. 
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FIG. 4: (Color online) As Fig. [Ubut for RHIC initial conditions. 
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FIG. 5: (Color online) As Fig. [2] but for RHIC initial conditions 






FIG. 7: (Color online) Evolution of scaled hadron densities n^/T 3 = YJs/T 3 ~ Y^ q s/T 3 (top panel) 
and related energy densities (bottom panel; the top curve is for the electro-weak contribution (7, 
i/'s, e, fj.) + pions (ir); the pion component is also depicted separately). Instead of p and p, the 
densities of nucleons (TV, n + ) and antinucleons (N, n_) are exhibited. For r\ = 10~ 10 . 
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FIG. 8: (Color online) Adiabatic path of cosmic matter in the phase diagram of strongly interacting 
matter for r\ = 10 -10 (solid curve for confined matter, h e ff = 10). The straight upper section 
(dashed curve) is for an approximation of deconfined matter. (The crossing does not necessarily 
imply phase equilibrium. For the given approximations, the turn from deconfined to confined 
adiabatic paths is by a Maxwell like construction with mixed phase; these path sections are not 
displayed.) The asterisks depict chemical freeze-out points from 
part, and LHC estimate mentioned in text there). 
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